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1  Introduction 

Markov  random  fields  based  on  the  lattice  Z 2  have  been  extensively  used  in 
image  analysis  in  a  Bayesian  framework  as  a-priori  models  for  the  intensity 
field  and  on  the  dual  lattice  (Z2)*  as  models  for  boundaries.  The  choice  of 
these  models  has  usually  been  based  on  algorithmic  considerations  in  order  to 
exploit  the  local  structure  inherent  in  Markov  fields.  No  fundamental  justifi¬ 
cation  has  been  offered  for  the  use  of  Markov  random  fields  (see,  for  example, 
GEMAN-GEMAN  [1984],  M ARRO QUIN-MITTER-P OG GIO  [1987]).  It  is 
well  known  that  there  is  a  one-one  correspondence  between  Markov  fields 
and  Gibbs  fields  on  a  lattice  and  the  Markov  Field  is  simulated  by  creat¬ 
ing  a  Markov  chain  whose  invariant  measure  is  precisely  the  Gibbs  measure. 
There  are  many  ways  to  perform  this  simulation  and  one  such  way  is  the 
celebrated  Metropolis  Algorithm.  This  is  also  the  basic  idea  behind  Stochas¬ 
tic  Quantization.  We  thus  see  that  if  the  use  of  Markov  Random  fields  in 
the  context  of  Image  Analysis  can  be  given  some  fundamental  justification 
then  there  is  a  remarkable  connection  between  Probabilistic  Image  Analysis, 
Statistical  Mechanics  and  Lattice-based  Euclidean  Quantum  Field  Theory. 

1This  research  has  been  supported  by  the  Air  Force  Office  of  Scientific  Research  under 
grant  AFOSR  89-0276  and  by  the  Army  Research  Office  under  grant  ARO  DAAL03-86- 
K-0171  (Center  for  Intelligent  Control  Systems) 
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We  may  thus  expect  ideas  of  Statistical  Mechanics  and  Euclidean  Quantum 
Field  Theory  to  have  a  bearing  on  Image  Analysis  and  in  the  other  direc¬ 
tion  we  may  hope  that  problems  in  image  analysis  (especially  problems  of 
inference  on  geometrical  structures)  may  have  some  influence  on  statistical 
physics. 

This  paper  deals  with  the  issues  just  described  and  above  all  it  suggests  a 
program  of  research  for  dealing  with  the  fundamental  issues  of  probabilistic 
image  analysis. 

What  is  the  fundamental  problem  of  Image  Analysis?  It  may  be  stated  as 
follows:  Given  noisy  images  (possibly  stereo)  of  the  visual  world  in  motion, 
represent  and  recognize  “structures”  in  some  invariant  manner.  For  example, 
if  the  structures  of  concern  are  three-dimensional  rigid  objects  then  we  want 
this  recognition  to  be  invariant  under  the  action  of  the  euclidean  group.  Thus 
symmetries  play  an  important  role  in  the  whole  process.  Ultimately,  we  want 
to  identify  the  mechanisms  (circuits)  which  perform  the  full  recognition  task. 
There  are  three  aspects  to  this  problem:  a)  the  identification  and  representa¬ 
tion  of  a  priori  knowledge  of  the  visual  world  without  which  the  recognition 
cannot  take  place;  b)  the  extraction  of  knowledge  from  the  imprecise  images 
available  about  the  visual  world  and  c)  understanding  the  correct  interaction 
between  a  priori  knowledge  and  the  imprecise  data  available  in  the  form  of 
images.  It  is  this  formulation  of  the  problematique  that  suggests  a  Bayesian 
framework.  In  this  paper  we  deal  only  with  some  partial  aspects  of  b)  and 
c)  and  it  is  here  that  the  connection  with  statistical  physics  enters. 

The  visual  world  is  a  continuous  world  and  if  we  believe  that  it  is  im¬ 
portant  to  capture  symmetries  then  mathematical  formulations  should  not 
be  on  the  lattice  C  but  should  be  based  on  R2.  Even  if  we  are  interested 
in  algorithms  on  a  digital  machine  which  necessarily  involve  discretization, 
then  in  order  to  capture  the  symmetries  in  the  limit  of  lattice  spacing  going 
to  zero  attention  must  be  paid  to  the  discrete  formulations  of  the  problem 
(see  for  example,  KULKARNI-MITTER-  RICH  ARSON  [1990]). 

This  paper  is  organized  as  follows.  In  section  2  we  discuss  Markov  fields 
and  Euclidean  fields  and  state  the  conjecture  regarding  Osterwalder-Schrader 
fields.  Section  3  is  concerned  with  a  variational  problem  in  Image  Analysis 
and  its  probabilistic  interpretation  and  it  shows  how  these  ideas  are  related 
to  those  of  Section  2.  Finally  in  Section  4  we  discuss  Stochastic  Quantization. 
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2  Markov  Fields  and  Euclidean  Fields.  (NEL¬ 
SON  [1973]) 

Let  Rd  denote  Euchidean  d-dimensional  space  ( d  >  2)  and  let  £(Rd)  denote 
Schwartz  space.  If  V  is  a  topological  vector  space,  a  linear  process  over  V  is 
a  stochastic  process  (p  indexed  by  V  which  is  linear  and  such  that  if  fa  —*  f 
in  V  then  <p(fa)  —*■  <p(f)  in  measure.  This  implies  that  <p(f)  for  /  e  V  is  a 
random  variable  over  (0,  E,  p). 

Let  ip  be  a  linear  process  over  <S(Rd).  If  A  C  Rd  is  open,  «4(A)  be  the 
<7-algebra  generated  by  <p(f)  with  support  of  /  C  A  and  if  A  is  any  subset  of 
Rd  let  -4.(A)  =  p]  *4(A')  where  A'  ranges  over  all  open  sets  containing  A. 

A'OA 

.4(A)  will  also  denote  the  set  of  all  random  variables  which  are  measurable 
with  respect  to  .4(A).  Let  Ac  denote  the  complement  of  A  and  d A  denote 
the  boundary  of  A. 

A  Markov  field  on  Rd  is  a  linear  process  <p  over  «S(Rd)  such  that  when¬ 
ever  A  is  open  in  Rd  and  a  is  an  integrable  random  variable  in  .4(A)  then 
£(a|>t(Ac))  =  EJ(a|^4(5A)).  Here  E  denotes  conditional  expectation. 

Let  0(d)  be  the  Euclidean  group  of  Rd.  By  a  representation  T  of  0(d) 
on  some  probability  space  (fi,  E,  p)  we  mean  a  homomorphism  rj  — ►  T(rj)  of 
0(d)  into  the  group  of  automorphisms  of  the  measure  algebra  and  this  group 
acts  in  a  natural  way  on  random  variables. 

A  Euclidean  field  is  a  Markov  field  over  S(Rd)  together  with  a  repre¬ 
sentation  T  of  0(d)  on  (0 ,E,p)  such  that  V/  G  <S(Rd)  and  rj  €  0(d) 

(Covariance)  T(rj)<p(f)  =  <p(f  o  t?-1),  (2.1) 

and  if  p  in  the  reflection  in  the  hyperplane  Rd_l 

T(p)a  =  a,  a  €  A(Rd~l).  (2.2) 

Consider  the  mapping 

:  Sn(fi,  •  •  •  ,  fn)  —  E(ifi(f[)  •  •  •  ^(fn)) 

and  assume  it  is  continuous.  By  the  Schwartz  Kernel  Theorem,  there  exists 
a  distribution  Sn  G  5'(Rdn)  such  that  S„(/i,  •••,/«)  =  Sn(fi  <g>  •  •  •  0  fn)- 
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Nelson  proves  that  the  sequence  of  distributions  Sn  satisfy  euclidean  in¬ 
variance,  symmetry  and  Osterwalder-Schrader  positivity,  namely,  if  A  is  the 
half-space  xd  >  0,  so  that  dA  is  the  hyperplane  Rd-1  and  Ac  is  the  half  space 
xd  <  0  then  E[(T(p)a)a]  >  0  where  a  €  *4(A)  and  T(p) a  e  v4(Ac). 

An  example  of  such  a  euclidean  field  is  obtained  as  follows.  Let 
be  the  real  Schwartz  space,  let  m  >  0  and  let  H  be  the  real  Hilbert  space  com¬ 
pletion  of  ^(R^)  with  respect  to  the  scaler  product  <  g,  (—A  +  m2)-1/  > 
where  A  is  the  Laplace  operator.  Let  <p  be  the  unit  Gaussian  on  H,  i.e.  ip  is 
a  real  Gaussian  process  indexed  by  H  with  mean  zero  and  covariance  given 
by  the  scalar  product  on  H.  Extend  <p  to  the  complexification  of  H  by  lin¬ 
earity.  Restricted  to  <S(Rd),  ip  is  a  linear  process  over  it.  Nelson  proves  that 
p  is  a  Euclidean  field,  that  is  it  is  Markov  and  satisfies  euclidean  invariance, 
symmetry  and  Osterwalder-Schrader  positivity.  Indeed  the  Markov  property, 
covariance  property  and  reflection  implies  Osterwalder-Schrader  positivity. 

Non-gaussian  random  fields  can  be  constructed  using  Multiplicative  func¬ 
tionals  or  measure  transformations.  Let  <p  be  a  Markov  field  over  <S(Rrf)  with 
the  underlying  probability  space  (0,  J7, //).  We  say,  that  a  random  variable 
P  is  multiplicative  if  for  every  open  cover  {A,}  of  Rd,  there  exists  strictly 

positive  fa  in  A(A,-)  with  fi  =  rift-  We  shall  see  later  how  we  construct 

* 

P(<p)%  fields  using  these  ideas. 

To  see  the  connections  with  Gibbs  density  let  us  proceed  formally.  Let 

=  ^((Vv?(x))2  +  m2(<p{x))2) 
and  consider  the  formal  expression 

exp(—  /  Ho(x)dx)  n  d<p{x). 

J  i€R2 


The  rigorous  interpretation  of  this  expression  is  as  a  Gaussian  measure  pc 
on  <S'(R2)  with  mean  zero  and  covariance  C  =  (—A  +  m2)-1.  This  measure 
corresponds  to  the  free  euclidean  field.  Non-gaussian  measures  are  obtained 
by  considering  formal  expressions 
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where  P  is  an  even  polynomial  and  A  is  a  coupling  constant.  The  above 
corresponds  to  the  canonical  Gibbs  density  for  transverse  vibrations  of  an 
elastic  membrane  subject  to  the  non-linear  restoring  force  F  —  —m2p(x)  — 
A P'(p(x)),  (after  integration  over  momentum  variables  p(x)),  and  '  denotes 
derivative). 

Now  if  p  is  the  unit  Gaussian  process  over  <S(R2),  then  P(p),  for  an  even 
polynomial  does  not  make  sense.  To  fix  ideas  let  us  consider  the  case  where 
P(p)  =  pA.  It  turns  out,  however  that 

:  ‘r’4  :  (d)  =  J  g(x)  ■  ¥>(*)  •  dx  for  g  6  L'(~]L°° 

has  a  well-defined  meaning  in  a  limiting  sense  in  £2(fl,  F,  fi),  where  :  p(x)4  : 
is  an  element  of  the  4th  Homogeneous  class.  In  a  more  concrete  way,  :  p3  := 
p3  —  3 E(p2)$  and  this  is  well-defined.  Now  if  we  consider 


exp (-fg(x)  :  p4(x)  :  dx) 
£(exp(— /3-(x)  :  p4(x)  :  dx) 


then  0  is  a  multiplicative  random  variable  and  we  are  able  to  construct  the 
non-gaussian  measure 


dp  = 


exp(—  /  :  p4  :  dx) 

Eex p(—  /  :  p4  :  dx)  ^ C 


This  measure  is  the  so-called  ((/?4)2-measure. 

The  proof  that  (<£>4)2-measure  defines  a  Markov  field  is  surprisingly  diffi¬ 
cult  and  has  been  accomplished  only  recently  (see  ALBEVERIO,  HOEGH- 
KROHN  and  ZEGARLINSKI  [1989]).  It  depends  on  the  complicated  theory 
of  local  specifications,  related  Gibbs  states  and  cluster  expansions.  On  the 
other  hand  the  property  of  Osterwalder-Schrader  positivity  is  much  easier  to 
verify.  With  a  view  to  proving  the  Global  Markov  property  and  for  problems 
in  Image  Analysis  we  advance  the  following  conjecture. 

We  define  an  Osterwalder-Schrader  field  (0-S)  on  R2  to  be  a  linear  pro¬ 
cess  p  over  <S(R2)  which  satisfies  the  Euclidean  covariance  property  (2.1)  and 
Osterwalder-Schrader  positivity.  For  example  an  Osterwalder-Schrader  field 
may  be  obtained  by  considering  a  function  of  a  Markov  field.  We  conjecture 
that  every  Osterwalder-Schrader  field  can  be  obtained  as  a  function  of  a 
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Markov  Field.  We  call  such  a  Markov  field  a  Hidden  Markov  Field  for  the 
0-S  field. 

Let  us  see  this  conjecture  in  the  familiar  context  of  stochastic  processes. 
Let  €  R-),  '■  D  — > ►  V  be  a  stochastic  process,  which  is  continuous  in 

probability,  stationary,  symmetric  (i.e.  tpt  and  are  stochastically  equiva¬ 
lent).  It  is  0-S  positive  if  V  0  <  t\  <  •  •  •  <  tn,  V  /  :  Vn  — ►  C,  bounded  Borel 
(V  is  a  topological  vector  space)  we  have 

<  f(<Pt i)  ■  •  •  •  *  S ‘f-tn)  >L2(n,r,P)>  0. 

Now,  let  <f>t  :  0  — >  V,  be  a  stationary,  symmetric  Markov  process  and  let 
tf;  :  V  — ►  V  be  bounded  and  Borel. 

Define  a  new  process  by 

¥»t(w)  =  V>(&(w))- 

<Pt  is  called  a  Markov  extension  of  <pt.  ipt  is  symmetric  and  stationary  and 
satisfies  0-S  postivity  but  is  not  necessarily  Markov.  The  conjecture  would 
be: 

Does  every  stationary,  0-S  process  have  a  Markov  extension? 

This  conjecture  in  discrete- time  is  not  true  (cf.  ARVESON  [1986])  but  in 
the  form  stated  previously  may  be  true. 

(It  is  clear  that  every  bounded  Borel  function  of  a  symmetric,  stationary 
Markov  process  is  0-S  positive). 


3  Probabilistic  View  of  Image  Segmentation 

We  think  of  a  noisy  image  as  a  function  g  :  Cl  — ►  R  where  0  C  R2  is  a 
bounded,  open  set.  We  assume  g  £  L°°(0).  A  variational  formulation  of  the 
image  segmentation  problem  due  to  Mumford  and  Shah  (cf.  MUMFORD- 
SHAH  [1989])  is  as  follows.  Approximate  g  by  a  function  /  and  a  closed  set 
r  CO  such  that  the  following  energy  function  is  minimized: 

E(f,T)  =  P  f  \f -g\2ix+  f  \Vffdx  +  aH'(T).  (3.1) 

Here  /  is  required  to  be  in  W/1,2(0\r),  Hl{T)  denotes  the  1-dimensional 
Hausdorff  measure  of  T  and  /?,  a  are  positive  constants.  We  would  like  to 
give  a  probabilistic  interpretation  of  E(f,T)  by  considering 
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exp(-  J  \Vf\2dx  -  aH1(T)) 

as  a  Gibbs  density  with  respect  to  a  suitable  reference  measure  which  is  to 
serve  as  a  prior  measure  on  (f,T)  and  exp(—  /  \g  —  f\2dx)  as  a  likelihood 
function  of  g  given  /.  The  choice  of  the  energy  functional  is  dictated  by  the 
requirement  that  /  should  approximate  g  in  the  L2-sense,  it  should  be  smooth 
away  from  the  boundaries  T  and  the  total  length  of  the  boundary  should  be 
short.  Note  that  for  fixed  T,  exp(—  fn/r  |V f\2dx)  rLceO\r  df(x)  would  have 
a  rigorous  interpretation  as  a  Gaussian  measure  pc  with  mean  zero  and 
covariance  C  =  (—A)  1  on  <S'(Q/r).  Using  the  recent  work  of  Surgailis  (cf. 
ARAK-SURGAILIS  [1989])  we  can  give  a  rigorous  interpretation  as  a  density 
to  the  following: 

exp(—  /  n  |A/|2cfo -2^£(rt) -y?(n)) 

jq\  U  r* 

t=l 

n 

which,  for  an  appropriate  <^(n),  is  a  density  over  /  £  Wq’2(Q\  [J  T,),  n  and 

t=i 

a  set  of  straight  lines  I\  which  form  together  with  dfl  a  polygonal  partition 
of  0  (cf.  MITT ER-ZEITO UNI  [1990]).  We  outline  here  the  basic  result  of 
Surgailis  which  constructs  a  measure  on  closed  polygonal  partitions  of  0. 

Let  0  be  a  closed  convex  subset  of  R2  with  smooth  boundary.  In  R2, 
choose  coordinates  ( t,x )  such  that,  for  all  y  £  fi,  t(y),  x(y)  >  0.  Let  Cq 
denote  the  lines  which  intersect  f2,  each  line  £  £  Cu  is  parameterized  by  its 
distance  from  the  origin  pi  and  the  angle  it  forms  with  the  t  =  0  axis,  at. 

Let  p(dt)  be  a  uniform  measure  on  the  set  a^pt  \i  G  Cq.  The  Poisson 
point  process  with  intensity  p{di)  will  be  denoted  and  the  measure  it 
induces  on  the  boundary  dQ  by  the  hitting  points  {( xt ,  ti)\ti  =  inf {t\£  gH}} 
is  again  a  Poisson  point  process  on  the  triple  (x,  t,  v)  with  intensity  gdQ.  Here 
and  in  the  sequel,  u  denotes  the  velocity  of  the  particle,  i.e.  the  tangent  of 
the  angle  formed  by  the  trajectory  and  the  t  axis. 

For  any  line  £  £  £n,  let  vt  denote  the  slope  of  £.  Clearly,  p(d£)  can  be 
considered  as  a  measure  on  vt  and  xe,  the  intersection  of  £  with  the  x  axis. 
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In  the  sequel,  we  consider  the  measure  p(dv,  dx )  obtained  from  the  uniform 
measure  p(d£)  on  a,  p,  i.e. 


p(d£)  =  p(dv,dx)  =  dx 


dv 


(1  +  i/2)* 


(3.2) 


Inside  Cl,  construct  a  point  process  on  the  quadruple  (t,  y,  v',  v")  with 
intensity 


pp(dt,  dy,  dv' ,  dv")  =  \v'  —  v"\dydt 


dv' 


dv" 


(1  +  ( l/')2)2  (1  -f  ( V")2)2 


(3-3) 


Finally,  construct  a  random  partition  of  Cl  as  follows: 

Pick  up  on  dCl,  n0  triples  ( t ,  x,  v)  according  to  the  law  pPu,  and  inside  Cl, 
n i  quadruples  ( t,y ,  v' ,  v")  according  to  the  Poisson  process  with  intensity  pp. 
At  each  of  those  points,  start  a  line  of  slope  v  (two  lines  of  slopes  v' ,  v"  in 
the  case  of  interior  points)  and  evolve  v  according  to  the  Markov  transition 
law 


P(vt+dt  £  du\vt  =  v)  =  |u  -  v 


dudt 

(l+«2)f ' 


(3.4) 


Finally,  at  each  intersection  of  lines  (when  viewing  it  in  the  direction  of  grow¬ 
ing  t )  kill  the  intersected  lines.  Clearly,  such  dynamics  describe  a  random 
partition  of  Cl  by  polygons,  c.f.  fig.  1.  The  basic  result  of  Arak  and  Surgailis 
is: 


Figure  1:  Random  Polygonal  Parition 
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Lemma  3.1 


P(n,£  £  d£i ,  *  =  1,  •  •  •  ,n)  =  —g(d£ i)  •  •  •  g(d£n)  exp(—  2  ^  C{£i))  (3.5) 

Til 


1  =  1 


where  £(£,)  denotes  the  length  of  the  i-th  segment. 

Note  that  due  to  the  presence  of  n  in  (3.5),  one  can’t  consider  (3.5) 
directly  as  a  candidate  for  a  density:  indeed,  if  one  were  to  consider  P(n,£,  £ 
£i  ±  e),  the  required  normalization  constant  (as  e  — »  0)  would  have  depended 
on  n  and  therefore,  a  path  with  no  jumps  will  be  infinitely  more  likely  than 
a  path  with  one  jump. 

One  way  out  of  this  problem  is  by  using  an  appropriate  definition:  Let 


Zn  =  f  ■■■(  P(n,£i  £  d£{,i  =  1,  •  •  •  ,n) 

JcD  JcD 

A  00 

z  =  £2”- 

n=0 

(4r)  is  the  probability  of  having  n  lines  in  a  specific  partition.  Now,  one 
may  define: 


Definition  The  prior  density  of  a  partition  (n,^t)  is  given  by 


/  n  \  rZnw  P(n,£i  t£i±e,i  =  1, 
p(n,({)  =  (T)hm - 5^ - 


,n) 


(3.6) 


Combining  these  ideas  with  some  of  the  ideas  contained  in  Dembo-Zeitouni 
(cf.  DEMBO-  ZEITOUNI  [submitted])  the  desired  result  referred  to  before 
can  be  obtained.  The  more  general  problem  of  interpreting  exp(—  /  |  V/|2da;— 
LT^F))  as  a  density  remains  open. 

The  interpretation  of  /? /n  \f  —  g\2dx  +  fn\r  |V/|2da:  +  aJT1(r)  as  a  pos¬ 
terior  density  accomplishes  something  important.  It  frees  us  from  obtaining 
Maximum  A  Posteriori  estimates  of  (/,  T)  via  minimization  of  the  above 
functional.  We  can  in  principle  obtain  other  estimates  such  as  conditional 
mean  estimates.  Indeed  for  closed,  convex  partitions  of  0  it  opens  up  the 
possibility  of  doing  inference  on  geometries  via  Monte-Carlo  simulations. 
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4  Stochastic  Quantization  (see  BORKAR- 
CH ARI- MITTER  [1988]  for  Stochastic  Quan¬ 
tization  of  (<£>4) 2  fields). 

The  problem  of  stochastic  quantization  for  the  energy  functional  (3.1)  is  to 
create  a  Markov  process  whose  invariant  measure  is  exp (—E(f,T)).  In  this 
section  we  wish  to  suggest  a  program  for  achieving  this  goal.  We  concentrate 
on  the  prior  density  exp(-  fn  jVfj2dx  -  aHx{ T)).  The  first  step  in  the  pro¬ 
cedure  is  to  replace  the  above  energy  functional  appearing  in  the  density  by 
an  approximating  functional 

/  \Vf\2(l-v2)ndx  +  a  [  ((\Vf\2  +  \Vv\2)(l  -v2)n  +  ^)dx. 

JQ  Ig 

(4.1) 

In  (4.1)  the  variable  u(x)  e  [0, 1]  should  be  thought  of  as  a  control  variable, 
which  controls  the  gradient  of  /  and  depends  on  the  discontinuity  set  T,  n  is 
a  parameter  which  tends  to  infinity.  The  above  expression  makes  sense  for 
fi  v  €  Wh2(Q,).  1  -(l-v2)n  approximates  smoothed  neighbourhoods  of  T  and 
the  approximate  boundaries  can  be  identified  with  (1  —  v2)n  ~  0.  If  we  denote 
by  Pn{B)  =  2(n  +  1)  Jb  y(l  —  v2)n\Vv\dx,  then  essentially  fin(Ct)  — ►  H4(T) 
in  a  weak  sense.  For  details  of  this  approximation  scheme  see  AMBROSIO- 
TORTORELLI  [1990]. 

The  first  step  in  the  program  would  be  to  identify 

exp (~£n(f,v)  JJ  d(f(x),v(x))  =  dpn 
x&K2 

as  a  measure  in  a  suitable  distribution  space.  The  covariances  involved  will 
have  appropriate  boundary  conditions.  One  can  then  study  the  weak  con¬ 
vergence  of  this  measure  d/xn  as  n  — *  oo.  The  functional  derivatives  of  Sn 
with  respect  to  /,  v  can  be  computed: 

jf  =  -V  ■  (V/  •  (1  -  u2)“)  (4.2) 

6-jk  =  — aV  •  (Vr  •  (1  —  u2)")  +  n(|V/|2  +  a| Vu|2)(l  —  v2)nv 
an 2 

+  Tg-o  (4.3) 
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(There  are  additional  terms  involving  the  normal  derivatives). 

In  analogy  with  the  study  of  stochastic  quantization  of  (<£>4)2  fields,  the 
problem  of  stochastic  quantization  for  fixed  n  is  the  study  of  the  coupled  pair 
of  stochastic  differential  equations 


dm  = 

S£n 

Sf 

■  dt  -f  dw(t) 

(4.4) 

dv(t)  — 

S£n 

8v 

•  dt  +  dv(t ) 

(4.5) 

where  w(-)  and  u(-)  are  infinite-dimensional  independent  Brownian  motions. 

5  Conclusions 

The  conceptual  program  outlined  in  this  paper  is  quite  general  and  may  be 
applied  to  other  variational  problems  arising  in  Image  Analysis,  for  exam¬ 
ple,  those  involving  curvature  terms.  These  variational  problems  may  be 
important  to  obtain  non-overlapping  segmentations  of  images  with  a  view  to 
identifying  occluded  regions. 
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